function E = ComputeStrain(F)
N = size(F, 1);
E = zeros(N, 2, 2);

for ii = 1:N
    F_ii = squeeze(F(ii,:,:));
    E(ii,:,:) = 1/2 * (F_ii.' * F_ii - eye(2));
end

end

